Quantum transport and spin dynamics on shearless tori 
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We investigate quantum dynamics in phase-space regions containing "shearless tori". We show 
that the properties of these peculiar classical phase-space structures — important to the dynamics 
of tokamaks — may be exploited for quantum information applications. In particular we show that 
shearless tori permit the non-dispersive transmission of localized wavepackets. The quantum many- 
body Hamiltonian of a Heisenberg ferromagnetic spin chain, subjected to an oscillating magnetic 
field, can be reduced to a classical one-body "image" dynamical system which is the well- studied 
Harper map. The Harper map belongs to a class of Hamiltonian systems (non- twist maps) which 
contain shearless tori. We show that a variant with sinusoidal time driving "driven Harper model" 
produces shearless tori which are especially suitable for quantum state transfer. The behavior of 
the concurrence is investigated as an example. 
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For a Hamiltonian system, the onset of classical chaotic 
dynamics is associated with the disappearance of phase- 
space barriers termed invariant tori: when the last invari- 
ant torus disappears, the chaotic diffusive motion is un- 
bounded ("global"). For a wide class of classical dynami- 
cal systems, described by so-called "twist-maps" and ex- 
emplified by the all-important Standard Map there 
is a single threshold for this process, i.e., once the last 
invariant torus breaks, the dynamics is globally diffusive 
for all parameters above the threshold. However, for an- 
other class of dynamical systems described by "non-twist 
maps" , this is not generally the case. A new class of tori, 
termed "shearless tori" [2] can be present. The classical 
properties of these shearless tori are attracting consider- 
able attention, due in part to their possible relevance to 
improved confinement of fusion plasmas in tokamaks 0], 
but also because the properties of non-twist maps are less 
well studied. New dynamical phenomena such as separa- 
trix reconnection lead to the intermittent re-appearance 
of the shearless tori which, when present, separate differ- 
ent regions of the chaotic phase space. 

In a quite different context, studies of quantum dynam- 
ics of spin chains, for example of quantum state transfer 
and entanglement generation, now play a central role in 
the burgeoning field of quantum information [4[. There 
is also growing interest in potential applications of non- 
linear dynamics in quantum information. In particular, 
the realization that a many-body Hamiltonian can be 
analyzed with the dynamics of "one-body image" quan- 
tum and classical Hamiltonians [5|, [6[ has proved very 
insightful. For example, the Heisenberg ferromagnetic 
spin chain, in a pulsed parabolic magnetic field, has the 
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Standard Map as its classical image [6] , provided position 
coordinates in the many-body Hamiltonian are mapped 
onto momenta of the image system and vice-versa. If, 
instead, a pulsed sinusoidal external magnetic field is ap- 
plied and the canonical coordinates are not interchanged, 
the spin system maps rather onto one of the best-known 
non-twist maps: the kicked Harper map. The quantum 
dynamics of the kicked Harper map have already formed 
the subject of several studies @,@]- From the viewpoint of 
entanglement screening, recently, calculation of the con- 
currence showed that nonlinear resonances in a classically 
mixed phase space induce robust entangled states [9j]. 

However, to date, the quantum dynamics near shear- 
less tori, and their potential for quantum information 
applications have never been investigated. In the present 
study, we have calculated quantum transport along 
shearless tori and found it to be quite different from the 
dynamics along ordinary tori. In particular, we find that 
shearless tori provide non-dispersive transport, while or- 
dinary tori are dispersive. In other words, for a spin 
chain, a spin wavepacket initially localized with start- 
ing conditions close to an ordinary torus, will rapidly 
become completely delocalized along the whole length 
of the spin chain. In contrast, for starting conditions 
near a shearless torus, an initially localized Gaussian 
wavepacket largely preserves its shape for long times, but 
moves along the spin chain at a near constant rate. In the 
chaotic regime, the shearless tori provide isolated regu- 
lar 'channels' through the chaotic sea which can be used 
to propagate coherent spin states (approximately) non- 
dispersively. We show that, for a driven Harper model, 
one may even start the wavepacket with zero momentum. 
We attribute the non-spreading of quantum wavepack- 
ets near shearless tori to a local quasienergy spectrum 
(i.e., for the Floquet states localized on the shearless 
tori) which is a perturbed harmonic spectrum. The corre- 
sponding evolution of the concurrence was also obtained. 

We consider the Heisenberg ferromagnetic spin chain 
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under an oscillating magnetic field. Here, we concentrate 
on a spin-1/2 and one-down-spin case, where a state is in 
a superposition of the states where one spin is down and 
the others are up. In that case, the model corresponds 
to a tight-binding model with an oscillating field of am- 
plitude Bo by the Jordan- Wigner transformation. Since 
the total S z is conserved, i.e. the one-down-spin state 
is conserved, the correlation terms disappear. Then, the 
Hamiltonian is written as 
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Here, ct and cj are the creation and annihilation opera- 
tors of a fermion at the jth site, respectively. For ferro- 
magnetic cases, J < 0. The total number of sites is N, 
and periodic boundary conditions are imposed. F(t) is a 
periodic function of time with period T. 

We solve spin dynamics from Eq. (pQ). However, to an- 
alyze the quantum behavior, we first note this: In the ab- 
sence of the time-dependent external field, the model (pQ) 
is associated with a dispersion relation, obtained from the 
Bethe ansatz, Ei = J(l — cos ki) with i — 0, 1, 2, . . . , N 
and where the E{ are energy eigenvalues and the k{ quan- 
tized wavenumbers. Equating the k{ with the momenta 
of a classical system [5|, |6|, i.e., k{ — > p and taking a con- 
tinuous position variable Xj = jj-j —> jf-x, we can see 
that Eq. (pQ) maps onto the one-body classical "image" 
Hamiltonian: 

H(x,p,t) = Jcosp + B F(t)cos[(27r/N)x]. (2) 

If F(t) = J2 n S(t — nT), we have the usual kicked Harper 
model. Here, we take F(t) = sinc^t, with uo = 27r/T, and 
investigate a modification we term the "driven-Harper 
model" . The corresponding classical equations of motion 
are given by: 



x = — Jsinp 

p= (2tt/N)B sin uot sm[(27r/N)x] 



(0<x< AT), 

( — 7T < p < 7r). 

(3) 

Integrating Eq. ([3]) numerically, we obtain x(t) and p(t). 
Classical surfaces of section (SOS) are obtained by plot- 
ting values of (x,p) at integer multiples of the period, 
t = n{2ir/uo), where n is an integer. Three SOS are 
shown in Fig.[Q Here, we set J = -1, B = 2, N = 100, 
and these values are fixed below. With decreasing fre- 
quency uj = 0.20, 0.16, 0.12, the classical phase space 
becomes increasingly chaotic; at uo = 0.16 all tori have 
disappeared; however, at uj = 0.12, shearless tori have 
re- appeared, embedded in two regular 'channels' running 
through a fully chaotic phase space (bar a few very small 
stable islands). 

For the kicked Harper map, shearless tori occur for 
p w ±7r/2 [2]. They correspond to an extremum of the 
rotation number; here || = at p ±7r. The shearless 
tori for uo = 0.12 are extremely distorted and span the 
full range of the map, from p = to ±7r/2. The kicked 
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FIG. 1: (Color online) Classical phase space for (a) uo — 0.20 
(near integrable), (b) uo — 0.16 mixed phase-space with no 
apparent tori, and (c) uo — 0.12 showing the re-appearance of 
shearless tori: phase space is almost entirely chaotic but reg- 
ular channels, containing shearless tori, which separate dif- 
ferent chaotic regions, are seen. The colored patched areas 
correspond to the initial wave packet of quantum dynamics 
in Fig. [2] 



Harper map, on the other hand, does not produce shear- 
less tori which include the p = region, so for this reason, 
using a driven rather than a kicked spin chain is advanta- 
geous. In other words, to prepare a spin- wave packet on 
the shearless torus region, one should excite low energy 
E(kj) c± spin- waves (i.e., kj c± since this corresponds 
to the bottom of the E{kj) energy band) in the region 
of the spin chain containing a node of the magnetic field 
(i.e., Xj w 7r/2 or Xj « 3ty/2). Exciting low-energy spin 
waves is a simpler procedure than selecting a more pre- 
cise range of higher-energies, which would be required if 
the shearless tori are only found at p w tt/2 as is the case 
for the kicked system. 

Returning now to our spin-Hamiltonian, Eq. (pQ), we 
compare the time evolution of a spin wavepacket with 
initial conditions corresponding, in the image classical 
system, to shearless tori with those of normal tori. We 
consider the quantum spin distribution: 



p(j,t) = \(j\m)\ 2 - 



(4) 



where \j) is the state where the jth spin is down, and 
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FIG. 2: (Color online) Time evolution of a quantum spin- 
wavepacket, initially centered at jo = 25, showing dispersive 
propagation on normal tori and non-dispersive propagation 
on shear less tori, (a) oo = 0.20 and ko — 1.0; wavepacket pre- 
pared initially on shearless torus of Fig. H^a), but in the near 
integrable regime, (b) uj — 0.20 and ko — 0.0; initial state is 
on an ordinary torus of Fig. [TJa) and rapidly delocalizes along 
the length of the spin-chain, (c) uj — 0.12 and ko = 0.0; the 
spin wavepacket initially on shearless torus embedded in the 
chaotic sea of Fig. H^c) propagates largely without dispersion, 
but the effect is less pronounced than in the integrable case, 
(d), (e), and (f) correspond to (a), (b), and (c), respectively. 
They exhibit the form of wavepacket s after t = 0T, IT, and 
10T. 



\ip(t)) is calculated from Eq. (pQ). Our initial state \ip(t = 
0)) is a Gaussian spin wavepacket with the phase space 
area corresponding to the classical plot in Fig. [TJ The 
center of the wavepacket is denoted by jo and ko in the 
real and momentum space, respectively. The width of 
the wavepacket is set to be Aj = 5. 

In Fig. [21(a), we show the behavior (for uo = 0.20) of 
a wavepacket initially centered on jo — 25 and ko = 1.0 
(i.e., mapping onto a regular phase-space region in the 
image classical system including a shearless torus), the 
distribution shows very little derealization for as long as 
the numerics were run (20 periods). Despite the same 
frequency, the distribution in Fig. [2Kb), which initially 
centered on the region corresponding to normal tori in 
the classical map, rapidly delocalizes and there is no dis- 
cernible localization at t — 10T. However, Fig. [2](c) 
(for uj = 0.12) shows again little derealization of the 
wavepacket initially centered on jo = 25 and ko = 0.0 
(i.e., on a shearless torus). Numerical experiments show 
that an equivalent classical ensemble of particles around 
a shearless torus spread slowly: the lack of shear also 
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FIG. 3: Quasi-energy spectrum local to shearless tori, showing 
that it corresponds to a near harmonic spectrum in the near- 
integrable limit and a perturbed harmonic spectrum in the 
chaotic limit. The spectrum is plotted for \Am\ 2 > 10 -3 . The 
low resolution spectrum is a Gaussian-smoothed one. The 
gray broad bars are just an eye guide, and their interval is 
2tt/5. (a) Quasienergy spectrum for uj — 1.0 showing equal 
quasi-energy spacing, (b) Low resolution (smoothed) spec- 
trum near shearless tori embedded in chaotic sea for uj = 0.12 
showing the spectrum approximates a perturbed harmonic 
spectrum, (c) Low resolution spectrum near ordinary torus 
at uj = 0.20. 



produces a strong classical effect. 

One can attempt to understand the quantum behavior 
from its quasi-energy spectrum: as this is a temporally 
periodic system, it is the eigenstates of the one-period 
unitary evolution operator (Floquet states) m (x,t) and 
associated eigenphases (quasi-energies) e m which repre- 
sent the stationary states of the system. One should first 
emphasize that the non-dispersive behavior is not sim- 
ply a wavepacket revival effect [Toj| where a wavepacket 
may partially or totally regain its shape after some re- 
vival time, dependent on the underlying frequencies. In 
our case, the wavepacket does not ever lose its essential 
shape. The behavior is rather analogous to motion of 
a wavepacket on a harmonic potential: the wavepacket 
oscillates in position without spreading; if the shape of 
the underlying potential is different (so our Aj is not the 
width of its ground state) the wavepacket exhibits only 
a certain "breathing" in its width, demonstrated in [1Q| . 
The modest variability in the wavepacket shape seen in 
Fig. [2] would be compatible with such breathing. 

We have calculated the local Floquet spectrum for the 
initial Gaussian wavepacket, i.e., the overlap between the 
wavepacket and the Floquet states, A m = (0 m (x,t = 
0)|^o). The probability |A m | 2 is plotted in Fig. [3] as a 
function of quasi-energy. For the wavepacket on a shear- 
less torus in a near-integrable regime, shown in Fig.[3^a), 
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FIG. 4: Time evolution of the concurrence, contrasting entan- 
glement transport by wavepackets moving along normal tori 
(left panels) with the non-dispersive shearless tori (right pan- 
els), (a) uo — 0.20 showing broad small peaks corresponding 
to Fig. [2(b), and (b) uo = 0.12 showing periodic sharp peaks 
corresponding to Fig. [2(c) . 



a 'harmonic' ladder of equally-spaced eigenstates is evi- 
dent. In Fig. 13(b), the shearless torus embedded in the 
chaotic regime, the low-resolution smoothed spectrum 
shows a set of equally spaced peaks. The interval of the 
peaks is 27r/5. This corresponds to the fact that the wave 
packet returns almost the same position as the initial one 
after 5 cycles (t = 5T) . We propose that this represents a 
perturbed harmonic spectrum found for eigenstates local 
to shearless tori, in the near-integrable regime. In the 
chaotic regime [Fig. [3(b)], the eigenstates are no longer 
completely localized on the shearless torus region and 
overlap with the chaotic regions; conversely, eigenstates 



with most support from the chaotic regions overlap with 
the torus region. Consequently, the states no longer form 
a ladder of individual states, but their amplitude is split 
redistributed into a cluster of nearby states. Neverthe- 
less, the presence of these regularly spaced clusters is 
sufficient to hinder dispersion of the wavepacket for a 
considerable time. For the ordinary torus, however, the 
spacing is irregular and even with smoothing, and the 
low-resolution spectrum has no equal spacing as shown 
in Fig. [3(c). 

The entanglement is also transported when the Gaus- 
sian packet travels. We employ the concurrence Cij, 
which is a measure of the bipartite entanglement of two 
sites i and j [11]. We concentrate on the cases corre- 
sponding to Figs. [2(b) and[2(c) as we are interested in the 
transport of quantum information. Figure [4] shows the 
time evolution of the concurrence, 625,26 ? 650,51 ? 675,76, 
and 6100,1 fc> r (a) uo = 0.20 and (b) uo = 0.12. Here, we 
should notice that the 100th and 1st sites are neighbor- 
ing sites because of the periodic boundary conditions. 
In Fig. [4(a), no sharp peaks appear, although we can 
see some small peaks. On the other hand, in Fig. [4(b), 
the concurrence has sharp periodic peaks. Moreover, 
several peaks are much higher than the initial value, 
625,26 (t = 0) — 0.22. This behavior of the concurrence 
corresponds closely to that seen in the spin dynamics. 

In conclusion, we find that quantum transport using 
shearless tori has a very different character from that as- 
sociated with normal tori. The quantum spin chain un- 
der an oscillating field can, for certain conditions, provide 
approximately non-dispersive transport even in a largely 
chaotic regime. The non-dispersion of Gaussian states is 
a property associated with large stable islands, in fact. 
These have been shown recently to preserve entanglement 
against noise [9]. The shearless tori may well provide a 
similar advantage while also transporting entanglement 
around the spin chain. 
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